Breathing oscillations of a trapped impurity in a Bose gas 
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Motivated by a recent experiment [J. Catani et at. larXiv:1106. 08281 /1 preprint, 2011], we study 
breathing oscillations in the width of a harmonically trapped impurity interacting with a separately 
trapped Bose gas. We provide an intuitive physical picture of such dynamics at zero temperature, 
using a time-dependent variational approach. In the Gross-Pitaevskii regime we obtain breathing 
oscillations whose amplitudes are suppressed by self trapping, due to interactions with the Bose 
gas. Introducing phonons in the Bose gas leads to the damping of breathing oscillations and non- 
Markovian dynamics of the width of the impurity, the degree of which can be engineered through 
controllable parameters. Our results reproduce the main features of the impurity dynamics observed 
by Catani et al. despite experimental thermal effects, and are supported by simulations of the system 
in the Gross-Pitaevskii regime. Moreover, we predict novel effects at lower temperatures due to self- 
trapping and the inhomogeneity of the trapped Bose gas. 



PACS numbers: 67.85.-d,03.75.-b,67.85.De 



I. INTRODUCTION 

The ability to trap and cool atoms of different species 
to ultra-low temperatures has led to the realisation of var- 
ious theoretical models in which the intriguing physics of 
binary mixtures of spin hyperfine states or different ele- 
ments can be studied [iHg. In particular, highly imbal- 
anccd mixtures have made it possible to investigate the 
dynamics, interactions and decohcrence of single atoms, 
generally referred to as impurities, immersed in a back- 
ground atomic gas [Mil- As a prominent example, sig- 
natures of polaron effects, caused by impurity-induced 
density fluctuations of the background gas, have been 
investigated theoretically [l^, [l!] and observed in exper- 
iments fillil. 



More recently, Catani et al. created a harmonically 
trapped impurity suspended in a separately trapped Bose 
gas They studied the dynamics of the system fol- 

lowing a sudden lowering of the trap frequency of the 
impurity. Primarily, they observed breathing oscillations 
of the width a of the impurity density distribution for 
various impurity-Bose gas interaction strengths. Several 
features of the experiment were amenable to interpreta- 
tion in terms of a quantum Langevin equation in conjunc- 
tion with a polaronic mass shift, however, we show that 
at even lower temperatures a different model is required 
to fully describe the dynamics of the system. 

In this article, we develop a versatile analytical model 
describing the experiment realised by Catani et al. at 
zero temperature. At first, our analysis is based on a vari- 
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ational approach in the Gross-Pitaevskii (CP) regime. 
We show that the impurity density distribution, which we 
describe by a Gaussian, has a width obeying a Newtonian 
equation of motion for a fictitious particle with position 
a. The potential governing this motion accounts for the 
quantum pressure of the impurity, the inhomogeneity of 
the trapped Bose gas and the localised distortion of this 
background gas induced by the impurity. The latter leads 
to a strong confinement of the impurity, known as self- 
trapping jl7H20| . Subsequently, we extend our model by 
including excitations of the Bose gas in the form of Bo- 
goliubov phonons. A variational ansatz, which describes 
the bosons as a product of coherent phonon states, allows 
us to track the evolution of the system, including the ex- 
change of energy between the impurity and the phonon 
bath. The dynamics turns out to be non-Markovian be- 
cause of the back-action of phonons created by the im- 
purity, and we show that the timescale of memory effects 
can be varied by adjusting the trapping parameters, al- 
lowing for a comprehensive study of the transition be- 
tween Markovian and non-Markovian dynamics. 

Our results reproduce the main features observed in 
the experiment of Catani et al. despite neglecting finite 
temperature effects, but also demonstrate that cooling 
the system further should introduce novel effects. We 
obtain an amplitude reduction of the breathing oscilla- 
tion for large impurity-Bose gas interactions. In particu- 
lar, we observe a sudden quench of the breathing ampli- 
tude once the impurity-Bose gas coupling substantially 
exceeds the boson-boson interaction strength. Similar to 
the experiment, we find that the breathing oscillations 
are damped due to the net dissipation of energy into the 
Bose gas. Moreover, our model explains the different be- 
haviour of attractive and repulsive impurities as a conse- 
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qucncG of the inhomogcneity of the trapped Bose gas. 



II. MODEL 

The specific system we consider consists of a single 
impurity of mass rua in a harmonic trap with frequency 
fla, and identical bosons of mass mt,, trapped by a sepa- 
rate external potential V{,(r), whose densities weakly in- 
teract with strength g. The impurity interacts with the 
Bose gas also via a density-density interaction of strength 
r]g, where the dimensionless parameter 77 controls the 
relative strengths of interactions. As such, the many- 
body Hamiltonian for the combined system is given by 
H = Ha + Hb + Hab, where 



Hb 



a = J drx^ 
drcp' 



''Kr X, 



2m,, 



Vb - 



Hab = V9 J drx'^X'^''^^, 



(la) 
(lb) 
(Ic) 



describe the impurity, bosons and impurity-Bose gas cou- 
pling, respectively. Here, x{r, t) and (/)(r, t) are the impu- 
rity and boson field operators. The single impurity has 
wavefunction x(r,i). 



III. BREATHING OSCILLATIONS IN THE GP 
REGIME 



(suppressed) from its decoupled value at the location of 
the impurity. The self-consistency of our assumption 
regarding the time-evolution of the condensate requires 
7io(r) 3> ??|x(i")P: ^-C-j sufficiently weak impurity-boson 
coupling. 

As a consequence of these approximations, Eq. (Pa|) be- 
comes the self-focussing non-linear Schrodinger equation 

We now focus on the important case of a spherically 
symmetric system in d dimensions with the Bose gas 
trapped in the harmonic potential Vb{r) = ^nibfllr'^ with 
frequency fib- The chemical potential in this case is 



Hb = fif^fc 



(2 + d)r(H-d/2) Ng 



2(27r)''/2 



hfibw^ 



2/(2+d) 



where Wb = ^h/rub^b is the width of the single particle 
ground state and V{p) = dxc~^x'P~^ is the Gamma 
function. 

We solve Eq. ^ within the variational ansatz 

X(r, cr, 7) = (ttct^) ''^'^ cxp {~r^ /2a'^ - i-ir"^) , (4) 

with time-dependent width a and phase 7. This ansatz 
was shown to describe a Bose gas in a harmonic trap ac- 
curately (see [2^ and references therein). We find equa- 
tions of motion for the Gaussian parameters by extremiz- 
ing the action S = J dt dr £ with the Lagrangian density 



In the GP approach, we replace (/)(r, t) by the conden- 
sate mode of the Bose gas Lp{Y^t). For low temperatures 
and in a dilute regime it follows from Eqs. ([T]) that the 
impurity and condensate wavefunctions evolve according 
to 
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respectively, where ip is normalised to the number of par- 
ticles N in the Bose gas f2l| . 

We consider the condensate in the Thomas-Fermi 
regime, in which the kinetic energy of the condensate, 
represented by the first term on the right hand side 
of Eq. ((2b)) . is negligible. If the impurity and Bose 
gas arc decoupled the condensate evolves as (po{r,t) = 
\/no(r)e~''' where the decoupled condensate density 
is given by no(r) = [fib - Vb{r)]/g for fib > Vb{r) and 
zero otherwise, and fib is the chemical potential. Return- 
ing to the case of non-zero impurity-boson coupling, we 
assume a similar form for the time-evolution of the con- 
densate f{r,t) — ■\/n(r)e~*'' It then follows from 
Eq. (|2bp that, for attractive (repulsive) interactions, the 
condensate density n(r) = no(r) — ?7|x(r)P is enhanced 



C = x* 



ih 



dt 2ma 



As a result, 7 ~ —niaO I2Tig and the width of the impu- 
rity (J obeys the equation of motion 



dVia) 

rriaa = 

da 

with the potential V ^ Vo + Vst + Vi„h, where 
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(6a) 
(6b) 

(6c) 



Here, R = y/2fib/mb^l is the Thomas-Fermi radius of 
the condensate and r{p,z) = [r{p)]~^ dxc~^xP~^ is 
the normalised lower incomplete Gamma function. 

The two terms in Vq describe the quantum pressure and 
the harmonic trapping of the impurity. The contribution 



3 




FIG. 1: One- dimensional homogeneous Boss gas. (a) The potential V governing the evolution of the width a of an impurity 
trapped with frequency 0,2 for various interactions strengths rj. The crosses mark the width a at the minimum of V for tight 
trapping with f2i = 5^2 ■ (b) The maximum a+ and minimum a- width during breathing oscillations as a function of 77. The 
analytical results of the GP approach predict that self-trapping strongly suppresses breathing oscillations (full blue) . Additional 
suppression occurs in the Bogoliubov approach due to dissipation into phonon modes (dashed orange) , matching the numerical 
results (green dots), (c) The frequency of the breathing oscillations as a function of 77. The parameters are set to mi, = 2ma, 
no = \/nih0.2lh and g = hO.2 ^/h/mb02. 



Vst is second order in rj and represents self-trapping ef- 
fects, which result from the deformation of the conden- 
sate due to the interaction with the impurity. The inho- 
mogeneity in the condensate density due to its trapping 
gives rise to the potential 14nh that is first order in 77. It 
has a simple form: The expression in the square brackets 
of Eq. ([5c)) is a monotonically decreasing function in a/R, 
taking the value 2/d at a/R = and asymptotically ap- 
proaching zero as cr/i? — !■ 00. Its gradient is largest for 
widths a ^ R. 

We see that Eq. ([5|), governing the evolution of the 
width (J, has just the form of Newton's second law for the 
motion of a fictitious particle with position a and mass 
nia in the potential V{a). This analogy allows us to un- 
derstand the behaviour of the width a as the impurity 
undergoes a procedure similar to that in the experiment 
by Catani et al. [l6j . In this procedure the impurity ini- 
tially feels a tight trap with frequency fii and the coupled 
impurity-Bose gas system is cooled to the ground state 
(although in the experiment such low temperatures were 
not achieved). At time t = Q the trap frequency is re- 
duced to and the subsequent evolution of the width 
for t > is observed. In our model this corresponds to 
the width of the impurity sitting at the minimum of V 
with frequency Via = ^1 for times i < 0. Then for t > 
the width moves in the potential V with lower trapping 
frequency ila = ^2- 

Our approach, with V = Vq, exactly captures the 
correct evolution of the impurity wavefunction for zero 
impurity-boson coupling. The square of the width oscil- 
lates harmonically with frequency 20,2 according to 



where ctj. = w'^{ni/il,2)^^ and Wa = \/hJm^Sh. is the 
width of the harmonic oscillator ground state. 



A. Homogeneous condensate 

For non-zero impurity-boson coupling, we first con- 
sider the special case of a homogeneous condensate whose 
decoupled density rip is uniform, corresponding to the 
regime of a flat trap i? ^ cr. In this case the potential 
14nh is also uniform and can be neglected, leaving the 
potential Y symmetric with respect to attractive and re- 
pulsive interactions. 

The self-trapping Vst imparts a force on the width 
of the impurity towards smaller values of a, shown in 
Fig. [TJa), and this has several effects. First, the mini- 
mum of the potential V is shifted to smaller widths. Sec- 
ond, it shortens the distance between the minimum of V 
for t < and its minimum for i > 0. Third, the curva- 
ture of potential V near its minimum is increased. These 
effects lead to an increased localisation of the impurity 
in the ground state at times t < 0. Further, at times 
t > the amplitudes of the breathing oscillations are 
reduced, along with an increase in the frequency of the 
oscillations. This is shown in Figs. [TJb) and[TJc), where 
wc have plotted the maximum and minimum widths, and 
the frequency of the breathing oscillations, respectively. 
The frequencies w are found from the best fitting curve 
of the form a'^{t)/w'^ = ^cos(a;i) -|- CTq to the oscillations, 
where A and cto are fitting parameters. 



B. Condensate in a harmonic trap 

In the general case the decoupled condensate density 
iM) depends on the radial position r and the potential V is 
augmented by Mnh- The total potential V resulting from 
a trapped condensate is plotted in Fig. [2ja). Since the 
contribution of Vlnh is first-order in 77 the symmetry be- 
tween attractive and repulsive impurity-boson coupling 
is broken. 
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FIG. 2: One- dimensional Bose gas in a harmonic trap, (a) The potential V governing the evolution of the width a of an 
impurity trapped with frequency 0.2 for positive (full) and negative (dotted) interactions strengths 77. The crosses mark the 
width o at the minimum of V for tight trapping with Q.i = 5f22- (b) The maximum (j+ and minimum ct_ width during 
breathing oscillations as a function of ri for a trapped (full red) and homogeneous (dashed blue) Bose gas. For a trapped Bose 
gas, breathing oscillations increase in amplitude with rj and at a critical value, here 77 « 3.1, undergo a sudden quench, (c) The 
frequency of the breathing oscillations cj as a function of rj. The parameters are set to ruf, = 2ma, TV — 200, Q,b ~ 0.2/^/^ and 
g = hQ,2 ^h/mhQ.2. 



For repulsive coupling, the force imparted by Vinh acts 
to increase the width of the impurity, opposing self- 
trapping. Specifically, in the regime Tj < 1 the inho- 
mogeneity of the Bose gas flattens the potential V near 
its minimum. In contrast, for sufficiently large 77 the self- 
trapping term dominates and the potential develops a 
narrow trough at a small width, as for the homogenous 
condensate. If interactions are attractive, Vlnh represents 
a force on the width of the impurity towards smaller val- 
ues of (J and enhances the self-trapping seen for the ho- 
mogeneous condensate. 

The effects of Vinh on the dynamics of the impurity fol- 
lowing a sudden decrease in Q.a are captured in Figs.[2Ib) 
and[2ljc). For weak repulsive interactions ?7 ^ 1, the ef- 
fect of the flattening of the potential is to increase the 
amplitude of oscillations and decrease the frequency. For 
strong repulsive interactions rj^ \ the creation of a sec- 
ond narrow minimum in the potential V results in small 
amplitude and high frequency oscillations. Our results 
predict a sudden transition between these two regimes, 
causing a sharp drop in the amplitude of oscillations and 
an increase in frequency. For attractive impurities, the 
effect of Vin\i is to enhance the effects seen for a homoge- 
neous condensate. 



IV. DAMPING AND MEMORY EFFECTS 

Up to now we have not included a mechanism for the 
transfer of energy between the impurity and the Bose gas. 
We now show, for a homogeneous Bose gas, that the os- 
ciUations of the impurity create Bogoliubov phonons and 
thereby lead to the exchange of energy with the Bose gas. 
Working in a volume V with periodic boundaries, we ex- 
pand = </7o+^0- The decoupled mode tpo only describes 
the Bose gas in the ground state if 77 = 0; therefore, in 
this picture, interactions with the impurity generate exci- 



tations even in the ground state of the system. Terms in 
H higher than second order in either the deformation Scf) 
or the coupling parameter 77 are neglected HI] . Then, 
expressing the deformation 5<j) = X]q["q('')^q + ^q('')^q] 
in terms of Bogoliubov modes 6q with momenta fi.q, the 
Hamiltonian of the Bose gas reads 

q 

Here, Eq is the energy of the mode tpo, Wq and 
fq solve the Bogoliubov-de Gennes equations, /iWq = 

■\/^q (^q + arc the phonon ener gies , and Cq ~ 

h?q^/2m}j are the free particle energies [2J]. Under the 
same approximations the interaction Hamiltonian simpli- 
fies to 

Hab = vgno + (^^q + ^q) /q ' C^) 

with the coefficients 

The first and the second terms in Eq. ([7]) represent the 
interaction of the impurity with the condensate mode ipo 
and the Bogoliubov phonons &q, respectively. Moreover, 
the latter describes the creation of Bogoliubov phonons 
due to a classical driving force /q. In our case, /q depends 
parametrically on the oscillating width a and thus the 
driving force is approximately periodic. 

We solve for both the ground state and evolution of 
the total system variationally within the ansatz I^P) = 
|cr, 7) 'Si \{c(q})- For the impurity, |(T, 7) corresponds 
to the Gaussian ansatz in Eq. The Bose gas is 

restricted to a product of coherent states |{aq}) = 
(g)^ g-lacil^/2gaqSij|Q^^ This is known to describe bosonic 
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FIG. 3: Damping of oscillations due to loss of energy, for weak 
coupling, (a) Damped oscillations of the width a of an im- 
purity in a trap with frequency 0,2 obtained analytically (full 
orange) and numerically (green dots). The trap frequency 
for t < is Sli = (b) The rate of energy exchange 

E between the impurity and the Bose gas. Positive (nega- 
tive) values correspond to energy gained (lost) by the Bose 
gas. The parameters are set to d = 1, 77 = 0.5, mt, = 2ma, 
no = Y^mbfiJ/ft and g — hQ.2^/hJriiiXh- The dashed line 
marks E = Q. 



modes coupled to a classical field such as a density and in 
particular Bogoliubov modes coupled to an impurity [l^l . 

We derive the equations of motion for the total system 
by extremizing the action S = J dt L with the Lagrangian 
L = {'^ \ (ihdt — H)\'^). Working in the thermodynamic 
limit, where V~^J2q ^ l/l^'"')''/ dq, we evaluate the 
integrals analytically using the small momentum approx- 
imation to the Bogoliubov dispersion relation cjq = cq, 
with the speed of sound given by c = a/ gno/mb- This 
approximation is accurate in the regime q ^ mi,c/h. We 
found no qualitative differences when proceeding numer- 
ically without making this approximation. 

Assuming the system is in the ground state at i < 0, 
the resulting equations of motion describe the width of 
the impurity a moving in the potential Vq augmented by 
a time-dependent force -Fphom representing the interac- 
tion of the impurity with the phonon bath, given by 



-^phon(^) 



S5(t,0) 



(8) 



c / dt 
Jo 



,c{t-t') [3I]2(f,t')_2c2(t-t')2]e 



1 /2 

with compound width T,{t,t') = [(T^{t) + cr^(i')] and 
constant K = 2dy/^/{i'!ry/'^r {1 + d/2)). The first and 
second terms in Eq. ([8]) correspond to the interaction of 
the impurity with phonons present in the initial state 
and those created during the evolution, respectively. 
These interactions introduce non-Markovian effects with 
a memory that decays exponentially on a timescale a/c. 
As a consequence, an experimentalist can, in principle, 
control the degree of non-Markovianity by altering g, uq 
and Wa- 

Due to the form of Fphon, the equation of motion for 
the width of the impurity a is an integro-differential equa- 
tion, which we solve using an iterative procedure: Ini- 
tially, we evaluate i^phon assuming a trajectory for a and 
integrate the equations of motion, subsequently we use 
the resulting trajectory to evaluate i^phon in the next it- 
eration. This procedure converges after a few iterations 
to a self-consistent evolution for a. The evolution of the 
width of the impurity after a sudden decrease in Qa is 
plotted in Fig. ^a). We find that the breathing oscilla- 
tions are damped out due to interactions with phonons. 
In Figs. Wljo) a-nd [TJc) we plot the maximum and min- 
imum widths, and frequencies as functions of coupling 
strength 77. Frequencies of oscillations uj are found from 
the best fitting curve a^(t)/wl = Ae^^^cos{u}t) + Ct + a'^ 
to our data, where A, B, C, A and ctq are fitting parame- 
ters. Dissipation before the first peak leads to a stronger 
suppression of maximum width than found when consid- 
ering self-trapping only. We find that dissipation also 
results in a greater dependence of frequency on coupling. 

To confirm our analytical approach we have numer- 
ically solved the coupled time-dependent GP Eqs. 
the details of which can be found in the appendix. For 
small I77I, we find good agreement on the frequency of 
oscillations and the width at the first maximum, while 
the analytics underestimates the heights of later peaks, 
as exemplified in Figs. [31[a). This leads to a close match 
of the numerically and analytically calculated maximum 
and minimum widths, and frequencies, as shown in Figs. 
[TJb) andfTfc). Also in these two figures, the numerics 
show that the symmetry between attractive and repulsive 
interactions is broken for large coupling; an attractive im- 
purity can infinitely enhance the condensate density at 
the origin, while an infinitely repulsive impurity only de- 
pletes the density of the Bose liquid at its location to zero 
(the Moses effect) [ll,[2^. 

We also calculate the energy deposited into the Bose 
gas E{t) = Eq^qdaqWI' - |aq(0)P). The rate of 
energy lost by the impurity is 

dE{t) AV.9'«o f cte-='*'/^'(*'°) 



dt 



S3(i,0) 



j:^t,t') 



/•* , t') - 2c2(t - t')2] e-'='(*-*')Vs^(^t') ] 

+Vo 1' 

which again is divided into two parts, representing in- 
teractions with phonons created at t < and t > 0, 
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respectively. In Fig. we plot this dissipation rate 
during the evolution shown in Fig. [3Ja). We find that 
most dissipation occurs when the width of the impurity 
is near its minimum value and that E can take negative 
values (the impurity can gain energy). It may, at first, 
appear at odds with the Landau criterion that dissipation 
is largest when the impurity density distribution is sta- 
tionary. However, such arguments do not apply here, as 
in [26j , since the impurity is not simply moving at a con- 
stant velocity. The form of E also differs markedly from 
that predicted by modelling the impurity as a damped 
harmonic oscillator, as done in [l6| . whereby dissipation 
is due to a friction force. The qualitative features of 
our analytic results are confirmed by numerical simula- 
tions, as shown in Fig.[3][b). Quantitative deviations arise 
from two approximations: We have neglected high order 
phonon effects and the impurity density does not always 
have a Gaussian form. 



and the impurity were 350 nK and 50 nK, respectively. 
We find that the typical energy scales of our system, 
hfli and the initial energy of the Bose gas, can be set to 
125 nK while retaining large oscillations rti/il,2 ^ 10 and 
the impurity width on the scale of micrometers. How- 
ever, mechanisms governing the evolution of the impurity 
at zero temperature, such as self-trapping, might still be 
relevant up to the temperatures realised in the experi- 
ment. 

Our flexible analytical framework allows for several 
natural extensions, e.g., impurities interacting directly 
and via the Bose gas, in the same trap or displaced. Fu- 
ture work will analyse many-impurity effects in such sys- 
tems. 
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V. DISCUSSION AND CONCLUSION 

Our results reproduce the main features of the experi- 
ment by Catani et al. [l6j even though we assumed zero 
temperature. Both the finite temperature experiment 
and our theory exhibit breathing oscillations of the im- 
purity after a sudden decrease in trapping frequency. In 
particular, for 77 = the oscillations observed in the ex- 
periment had a frequency and maximum to minimum 
width ratio close to the exact zero temperature values. 
Also, in both experiment and theory, large impurity-Bose 
gas coupling suppresses the amplitudes of breathing os- 
cillations and leads to damping. However, our zero tem- 
perature calculations predict a richer dependence of am- 
plitude on coupling, including a sudden quench, resulting 
from the interplay of self-trapping and inhomogeneity in- 
duced by the Bose gas trapping. We also predict a de- 
pendence of oscillation frequency on the impurity-Bose 
gas coupling, indicating that the lack of a dependence 
observed by Catani et al. is a finite temperature effect. 
Our intuitive model explains differences between repul- 
sive and attractive impurities, and highlights the way in 
which the non-Markovianity of the system arises and how 
it may be controlled. 

In general, we expect the zero temperature approxima- 
tion to be valid if the thermal energy of the initial state 
k^T is small compared to the relevant energy scales of the 
system at i < 0. In the experiment by Catani et al., prior 
to the equilibration the thermal energy of the Bose gas 
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Appendix A: Numerical simulations 

We solve the coupled time-dependent CP Eqs. ([2]) for 
a homogeneous Bose gas numerically using time-splitting 
methods. We assume for numerical ease that the conden- 
sate is enclosed in an infinite spherical box with a radius 
TZ that is large compared to the width a and the heal- 
ing length ^. More precisely, we choose TZ = AQwa and 
solve for tf{r) and x(r) on a spatial grid of 10^ points 
using timestep 10~'^r2^^. The ground states are calcu- 
lated using the normalised gradient fiow (27l - [29| and the 
time-evolved states are calculated by applying a method 
combining the time-splitting technique to decouple the 
nonlinearity [28l . |30| with the second order finite differ- 
ence and Crank-Nicolson methods to discretise the spa- 
tial and temporal derivatives, respectively 31[. To calcu- 



late the width of the impurity cr, we find the best fitting 
parameters to |x('')l = i?exp(— r^/2(T^). Maximum and 
minimum widths are calculated from this data. 
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